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Abstract 

We study the free energy of a particle in (arbitrary) high-dimensional Gaussian random potentials with 
isotropic increments. We prove a computable saddle-point variational representation in terms of a Parisi- 
type functional for the free energy in the infinite-dimensional limit. The proofs are based on the techniques 
developed in the course of the rigorous analysis of the Sherrington-Kirkpatrick model with vector spins. 



1 Introduction 

Recently, considerable (renewed) attention in the theoretical physics literature has been devoted to Gaussian 
random fields with isotropic increments viewed as random potentials, see, e.g, the works by Fyodorov 
and Sommers [8], Fyodorov and Bouchaud [7|, and references therein. In particular, it was heuristically 
argued in these works that Parisi's theory of hierarchical replica symmetry breaking (Parisi Ansatz, cf. 1 1 1 1) 
is applicable in this context. In the probabilistic context, these results provide rather sharp information 
about the extremes of the strongly correlated fields with high-dimensional correlation structures, which is 
a challenging area of probability theory fl4l l4l l2l [3l [T71 [181 . 

In this note, we initiate the rigorous derivation of the results of J8] [7) . We concentrate on the compu- 
tation of the free energy of a particle subjected to arbitrary high-dimensional Gaussian random potentials 
with isotropic increments. In the high-dimensional limit, we derive a computable saddle-point representa- 
tion for the free energy, which is similar to the Parisi formula for the Sherrington-Kirkpatrick (SK) model 
of a mean-field spin glass. Our proofs are based on the local comparison arguments for Gaussian fields 
with non-constant variance developed in [5 |, which are, in turn, based on the ideas of Guerra [9|, Guerra 
and Toninelli 1101 . Talagrand lfl6ll and Panchenko lfl3l . 

This note is organised as follows. We state our results in Section[2] The proofs are given in Sections [3] 
and|4] In Section|5] we give an outlook and announce some important consequences of the results of this 
note. In the Appendix, we provide some complementary information for the reader's convenience. 

2 Setup and main results 

Consider the Gaussian random field with isotropic increments X = Xn = {Xn(u) : u E R N }, N £ N. The 
adjective "isotropic" means here that the law of the increments of the field X is invariant under rigid 
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motions (= translations and rotations) in R . We are interested in the case N 1 and in the case of 
strongly correlated fields with high-dimensional correlation structure. Therefore, we assume that the field 
Xn satisfies 

E[(X N {u)-X N (v)) 2 ] =D^||m-v||^ =:D N (\\u-v\\l), u,v£R n , (2.1) 

where || -\\2 denotes the Euclidean norm on R N and the correlator D : R + — > R + is any admissible function. 
Complete characterisation of all correlators D that are admissible in < |2.1) , for all N, is known, see Theo- 
rem A. 1 Note that the law of the field Xn is determined by ( |2.1| > only up to an additive shift by a Gaussian 
random variable. In what follows, without loss of generality, we assume that Xn(0) = 0. 

We are interested in the asymptotic behaviour of the extremes of the random field Xn on the sequence of 
the particle state spaces Sn C R n as N t +°°- The state spaces are assumed to be equipped with a sequence 
of a priori reference measures {/%} C «^finite(>%)' We now define the main quantities of interest in this 
work. Consider the partition function 

Zrtifi) :=J s Md")exp(/3v^v(«)) . P G R. (2.2) 

We view ( |2.2[ ) as an exponential functional of the field Xn, which is parametrised by the inverse temperature 
[5. Heuristically, for large j3 (i.e., j3 t +°°), the maxima of the field Xn give substantial contribution to the 
integral ( |2.2[ ). The A^-scalings in ( |2.2[ i, ( |2.1[ ) and the "size" of Sn are tailored for studying the large-A^ limit 
of the quenched log-partition function: 

p N (P):=^logZ N (P), peR. (2.3) 

For comparison with the theoretical physics literature, let us note that there one conventionally substitutes 
/3 i — y —ft in ( |2.2| i (this has no effect on the distribution of Zn due to the symmetry of the centred Gaussian 
distribution of the field Xn), and considers instead of §2.3) the free energy 



fN(P):=-~P N (P), PeR+. (2.4) 

Assumptions. Informally, we require the particle state space Sn to have an exponentially growing in N 
volume (respectively, cardinality, if Sn is discrete). In particular, using physics parlance, this assures that 
the entropy competes with the energy (given by the random field Xn) on the same scale. More formally, we 
assume 

S N := S N , ScR. (2.5) 

Let ji £ ^finite(S) be such that the origin is contained in the interior of the convex hull of the support of ji. 
Define jj,N := /J G .^finiteOSV)- A canonical example is the discrete hypercube Sn '■= {— 1; 1} equipped 
with the uniform a priori measure, i.e., jx({u}) := 2~ N , for all u £ Sn- 

Parisi-type functional. To formulate our results on the limiting log-partition function, we need the 
following definitions. Given r £ R + , consider the space of the functional order parameters 

3£(r) :— {x : [Q;r] —>■ [0; 1] | xis non-decreasing cadlag, x(0) = Q,x(r) = 1}, (2.6) 

It is convenient to work with the space of the discrete order parameters 

SC^if) : = {x £ 3£{f) | x is piece-wise constant with at most n jumps}. (2.7) 

Let us denote the effective size of the particle state space by 

d := sup — sup ||w|| 2 . (2.8) 
n \N ueS N J 

For what follows, it is enough to assume that r £ [0;d] in ( |2.6| l. Note that, in case ( |2.5| >, d — sup M ^£ u . 

Now, let us define the non-linear functional that appears in the variational formula of our main result. 
We do it in three steps: 
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1. Given large enough M € R + , define the regularised derivative D'- M : R + — » R of the correlator Z) as 

D ^ (r): =< p>'M, ^[l/M; + oo), 
W |M, re[0;l/M). 

Given r,MfR + , define the function Or : r: r ► R as 

r (M) ( ? ):=^ M (2(r- ? )) + iD(2(r- ? )), ?6 [-r;r]. (2.10) 

2. Given r £ R+, i£ and the (regular enough) boundary condition h : R — > R, consider the semi- 
linear parabolic Parisi's terminal value problem: 

Uf(y,q) + \D'> M {2{r-q)) (d 2 q f(y,q) +x(q) (d y f(y,q)) 2 ) = 0, (y,q) G R x (0,r), ^ ^ 
\/(y,l)=%), yeR. 

Let f^l : [0; 1] x R + ->• R be the unique solution of (JZTTJ. Solubility of the Parisi terminal value 
problem ( |2.1 1[ ), its relation to the Hamilton-Jacobi-Bellman equations and stochastic control prob- 
lems is discussed in a more general multidimensional context in [5, Section 6]. 

3. Given the family of the (regular enough for ( |2.1 1| > to be solvable) boundary conditions 

:R->R| A eR}, (2.12) 
and given r 6 [0;d], define the local Parisi functional £?($,r,g): X (r) — > R as 

&(fi,r,g)\x] := lim ( inf \A (0, 0) - Xr] - %■ [' x(q)d0 ( r M) (q)) , x e %{r). (2.13) 



In ( |2.13[ ), the integral with respect to of 1 ^ is understood in the Lebesgue-Stiltjes sense. 

Main results. Let us start by recording the basic convergence result for the log-partition function. 

Theorem 2.1 (Existence of the limiting free energy). For any /3 > 0, the large N -limit of the log-partition 
function exists and is a.s. deterministic: 

Pn(P) > p(P)i almost surely and in l) . (2.14) 

In addition, for any iVeH the following concentration of measure inequality holds 

( Nt 2 \ 

P{\p N (P)-E\p N (P)]\>t}<2exp\-— -J, teR+. (2.15) 

The main result of this work is the following variational representation for the limiting log-partition 
function in terms of the Parisi functional ( |2.13| l. 



Theorem 2.2 (Free energy variational representation, comparison with cascades). Assume ( |2.5| >. Let the 

family of boundary conditions \2.\2\ be defined as 

gx{y) ■= log j s H(du)exp(puy + Xu 2 ), yeR. (2.16) 

Then, for all j3 G R, 

p(P) := sup inf (^(/3,r,,g)[x] — &(r)[x]) , almost surely and in L l , (2.17) 

r£ [Q-d]xsSC(r) 

where the remainder term S%{r) : 3£{f) — > R+ is afunctional on ^ if) taking non-negative values {see 
( |4.23| > for the definition). 
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The sign-definiteness of the remainder term &(r) immediately implies the following bound. 
Corollary 2.1 (Log-partition function upper bound). For all j3 G R, 

/?(J3) < SU P inf ^(P, r >g)[ x ]> almost surely. (2.18) 

Remark 2.1. In the case \A.A) , the field \2.2Q\ has a feature, which is not within the assumptions typically 
found in the literature ilPl 1701 176] 1751 \13\l : the correlator D is not of class C 1 , namely, D can have a singular 
derivative at 0. To deal with the singularity, we need a regularisation procedure, cf. \2.9\ and ( |2.13| l. 

Heuristics. It is natural to ask the following questions: Why is Parisi's theory of hierarchical replica 
symmetry breaking [ 1 1 1 (which is usually behind the functionals of the type ( |2.13| >) applicable to Gaussian 
fields with isotropic increments satisfying ( |2.1[ )? Where are the "interacting spins" in the present context? 

A hint is given by the following observation. Define 

1 N 

(u,v) N :=~Y, u i v i> «.veR w . (2.19) 

Let us fix r 6 [0; d] . By ( |A.6) , the restriction of the field Xn with isotropic increments to a sphere with 
radius r centred at the origin, leads to the mixed p-spin spherical SK model (cf. [ 15 1) with the following 
covariance structure 

E[X N (u)X N (v)]=D(r)--D(2(r- (u,v) N )) =: G r ({u,v) N ), \\u\\ 2 2 = ||v||f = rtf, (2.20) 
where G, : R+ — > R is given by 

G r {q):=D{r)-^D{2{r-q)), qeR+. (2.21) 

Thus, ( [230) impl ies that, given r, each field of the type ( |2.1[ ) induces a mixed p-spin spherical SK model 
with the convex correlation function G r (see Remark |A.2| >. It is this convexity that leads to the sign- 
definiteness of the remainder term in ( |4.24| i and allows for the proof (along the lines of lfT6ll ) of Theo- 
rems [5J] and for all admissible correlators. 

Our proof of Theorem |2.2| exploits the observation ( |2.20[ ) and combines it with the localisation tech- 
nique of 0. By means of the large deviations principle, this technique reduces the analysis of the full 
log-partition function ( |2.3[ ) to the local one, where ( |2.20[ > approximately holds true everywhere. The price 
to pay for this reduction is the saddle point variational principle ( |2.17| i, which involves the Lagrange mul- 
tipliers that enforce the localisation. 

3 Existence of the limiting free energy 

In this section, we prove Theorem |2.1| 

Proof of Theorem \2. 1 | Proof of p,15[ ). By Remark |A.l| we have 

Var[X N (u)]=D N {\\u\\l) <D{d), u G R N . (3.1) 

Therefore, the concentration of measure inequality ( |2.15| l follows from |5j Proposition 2.2]. 

Proof of the convergence ( |2.14| i. The result can be proved along the lines of 1 10 Theorem 1]. In [ 10, 
eq. (7)], it is assumed that the covariance structure of the random potential depends on the scalar product 
(overlap) of the particle configurations in a smooth way. Therefore, using the terminology of Remark |A. 1| 
only the short-range case is covered by [10, Theorem 1]. Indeed, in that case, the covariance of the field X^ 
satisfies ( |A.1| >, where the function B is analytic and convex, which follows from the representation ( |A.2| >. 
Therefore, [ 10, Theorem 1] is applicable with Qn{u, v) := TV -1 \\u — vH^, for u,v € R N . 

In the long-range case ( |A.6| >, the proof of the [ 1 1 requires some care, because the covariance structure 
of the field Xn (cf. ( |A.6| >) does not depend on the scalar product ( |2.19| > only, and, moreover, the correlator 
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D is not of class C 1 (cf. Remark 2.1 1. For the reader's convenience, we now retrace the main parts of this 
argument. Given E N, we prove the convergence of ( |2.14| i along the subsequences {A^ := N K }Ketti- 

Ak) 



Convergence along other subsequences then readily follows. Consider N independent copies {X^ 1 \kE 

(3.2) 



[N]} of the field X^ Kl ■ Given an interval V C [0;d], define the localised state space as 

S N (V) := {u e S N : e N-V}. 
Given a random field C = {Cn(u) \ u G R. N }, denote the corresponding local partition function by 



Z N (P,V)[C}:= [ fi N (du)exp(pVNC N (u) 

Js N (v) 



(3.3) 



In what follows, for u G R N , v G R M , we denote by u 1 1 v the vector in R N+M obtained by concatenation of 
u and v. Define the Gaussian field Y as 



y^( M ( 1 )n«( 2 'n...i li( M):=4=j;xg_ i («), M «GR^-S fcG[iV]. 
Due to independence, 



(3.4) 



Cov 



ii ii ... ii «W),y M (F(" i, V W ii ... ii v W ) 



= E Cov < , (» w ),4? , (^) . « ( V*> g r^- 1 , * g [iv]. 



Let us define 



Z Nk (P,V)[C] := / ^(dw)exp (pVNC N {u)) , 



where 



Sn k (V) := |m = ii « 



(2) 



ii...iim (a ' ) £S Nk : ||uW||^ GiVk-l-V, feG 



(*)ll? 



[N]}- 



Let us note that Sn k (V) C Sa^ (V), and, therefore, 

z NK (p,v)>z NK (p,v). 

The product structure (|3.7|i and independence (|3.4-|> imply 



(3.5) 

(3.6) 
(3.7) 

(3.8) 



Mr 



-E 



logZ NK (p,V)[Y N>K ] 



-E 



1 
1 



lognz^os.v)^] 

EflogZ^^,^)^]]. 



(3.9) 



For £ > 0, set V* := [ie; (/ + l)e], i G IN. By the Gaussian comparison formula |5 , Proposition 2.5], 



N K 



logZ NK (P,Vi)[X NK 
p 2 ri 



1 



E[\ogZ{p,Vi)[Y NJC l 



dt I % K {t){Au) I 95^(0 (dv) 

1 



VarX^ («) - - £ VarX^ , («« ) (3 



10) 



Cov[^( M ),^(v)] - - £ Cov [Xn^ (aW),^^*)) 



where ^v^(f) G .J#\(Sn k ) is the interpolating Gibbs measure with the density 

^^( M )=exp(j3^(V^(«) + \/r^7y^( M ))), K6%(Vi). (3.11) 
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Using ( |A.6| l, the smoothness of the correlator D on (0; +°°), the fact that D is non-decreasing, continuous 
at 0, and D(0) = 0, we get 



sup 

ues NK (Vi) 



iV fc=l 



<D(e), ieN. 



(3.12) 



As for the covariance terms, the concavity of the correlator D (cf., Remark A. 2 1 and the explicit covariance 
representation \A.6\ assure that 



sup ( Cov [X Nk (u),X Nk (v)} - - £ Cov X Nk _ { (vW) ) < D(e) 



Therefore, combining < |3l8] ), ([3T9J, (|37TOj, ( |3~T2] > and ( |3J3) we get 



(3.13) 



~ E [logZ^ (j8 , V,) pfrj] > - J-E [logZ^, 03 , V}) [X Nk J] - CD(e), ieN. (3. 14) 



The proof is finished by using the concentration inequality ( |2.15| > to remove the localisation in ( |3.14| i, as in 
ifTOl Theorem 1]. 

□ 



4 Comparison with cascades 



In this section, we prove Theorem 2.2 The proof follows the strategy that was previously implemented 
in UJ Section 5]. The appearance of the auxiliary structures below can be made more transparent by the 
"cavity" arguments, as is done in the seminal work of Aizenman et al. (T]. 



4.1 Auxiliary structures 

Consider the auxiliary index space si — si„ := N", n G IN. Let us define the projection operator si 3 a H> 
[<x]k '■= (05i, . . . , OCfc) G H k , for k G [n]. It is useful to treat the elements of si as the leaves of the tree of 
depth n. We use the convention that [a]o = 0, where denotes the root of the tree. Given a leaf a € srf ', 
we think of { [o^] '■ k E [n] } as of the sequence of branches connecting the leaf a to the root 0. We equip si 
with a random measure called Ruelle 's probability cascade (RPC). Let us briefly recall the construction of 
the RPC, see, e.g., [1J for more details. Note that each function x G can be represented as 



*(«) = L^ 1 [?i;?i+i)( r )' 

i=0 



where x — {xk}tto m< ^ 1 = fe}fc=o satisfy 



=: xq < x\ < . . . < x„ < x n +\ '■= 1, 
=: q < qi < ■ ■ ■ < q„ < q„+i := r. 



(4.1) 



(4.2) 



To define the RPC, we need only the sequence x as in ( |4.2[ ). Consider the family of the independent 
(inhomogeneous) Poisson point processes , | (X G si ,k G [n]} on R + with intensity 



3(4.t i r' r 'eK + , fee [l;n]nlN. 



(4.3) 



To each branch [tt]k, ex. G si, k G [n] of the tree we associate the position of the o^-th atom (e.g., according 
to the decreasing enumeration) of the Poisson point process ^k.[a] k _ { - The RPC is the point process RPC = 
RPC(jci , . . . ,x n ) := Y*ae.o/ <>RPC(a)> where RPC(a), a G si is obtained by multiplying the random weights 
attached to the branches along the path connecting the given leaf a G si with the root of the tree: 



RPC(a):=n&.[«]n(°*)- 

k=l 



(4.4) 
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Since L aG .0/RPC(o!) < °°, the RPC can be thought of as a finite random measure on stf with (abusing the 
notation) RPC({a}) := RPC(a), for a € si. To lighten the notation, we keep the dependence of the RPC 
on x implicit. 

Recall (3.2) , Given the sequence x as in ( |4.2| i and any suitable Gaussian field C := {C(m,o:) | u G 
5a?, 0£ € ,<z/}, let us define the extended log-partition functional <£>n(x,V) as 



^y)[C]:=iE 



log ( / n(du) [ RPC(da)expfj3v / iVC( M ,a) 

\Jsn(v) \ 



(4.5) 



where the RPC is induced by I. 

Let us use the remaining from the order parameter i£ bit of information, namely, the sequence 

q = {qk}'l=Q, as in (4.2) , to construct the Gaussian cavity fields indexed by Sn x si. To this end, define the 
lexicographic overlap between the configurations a' 1 ', a' 2 ' G si as 



Z(a (1) ,a (2) ):= 



0. 

max < 



(4.6) 



v ;$yk e[N]:[a {l \ = [a {2 ' , ] k ^ , otherwise. 

Let us define (slightly abusing the notation) the lexicographic overlap q : si 1 — ► [0; 1] as 

q{a {l \a^):=q l(aWaW) . (4.7) 

Given ^ as in (4.2) , the cavity field is the Gaussian field A = A^ — {An(u, a) \ u G Sjv, OC € .e/} such that 
Cov f 



A( M )(«,«W),AW(v,a^)] =D'' M (2(r-^(a( 1 ),a( 2 ))))( M ,v) A ,, a«,a( 2 W, 

(4.8) 

The existence of the cavity field A is guarantied by the following result. 



Lemma 4.1 (Existence of the cavity field). For any sequence q as in (4.2) and large enough M € R+, f/ier<? 
exists f/;e unique (in distribution) Gaussian field satisfying (4.8). 



Proof. Since the distribution of the Gaussian field is completely identified by the covariance, the unique- 
ness follows once we prove the existence. For this purpose, we first construct the Gaussian field a = 
{a( M )(a)Wwith 



Cov 



2 W (a (i) )jfl (M) (a (2) ) ] =D '^(2(r-q(a^\a^))) l a« a (2) e^. (4.9) 

To construct the field a^ explicitly, we define 

m k :=D lM {2{r-q k )), ke[n + l]. (4.10) 

The representations (A. 3) and (A.4) , guarantee that the sequence (4.10) is non-decreasing. Therefore, we 
can set 

a^(a) := £ (m k+l -m k f 2 gfX, 



(4.11) 



k=l 



where {g^ I a S & ' i k S W} are i.i.d. standard normal random variables. A straightforward check shows 
that the covariance structure of (4.1 1) satisfies (4.9) . 

To finish the construction, for i E [N], let a^= {a\ M \a)} ae ^ be the i.i.d. copies of the field a' M ' = 
W M \a)} ae ^. Define 



i N 

A { ™\u,a):= -=Y^ a i \ a ) u i^ ueS N , a est. 



An inspection shows that the field (4.12) satisfies (4.8) . 



(4.12) 

□ 
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4.2 Interpolation 



In this section, we shall apply Guerra's comparison scheme (cf. to the Gaussian field with isotropic 
increments satisfying \2.l) . To this end, we restrict the state space of a particle to a thin spherical layer. 
This assures that the variance of the field Xn does not change much. We refer to this procedure as local- 
isation. Then, we interpolate between the field of interest Xn and the cavity field ( |4.12[ ) and compare the 
corresponding local log-partition functions. We use the auxiliary structures from Section [4~T| 

Given x g ^'(r) and large enough M £ R + , let us consider the following interpolating field on the 
extended configuration space S^x si 



H( M \u,a):=VtX N (u) + VT^tA [ ™>(u,a), re[0;l], u e S N , 



AM) 



(4.13) 



where Ajy is the cavity field with ( |4.8| ). In the usual way, the field ( |4. 1 3] > induces the local log-partition 
function 



(p [ f\t,x,V):=<t> N (x,V)[H t }, Vc[0;d], xtSC'ij). 



At the end-points of the interpolation, we obtain 



(p^ } (O,x,V) = N (x,V)[A^] and (p ( N M) (l,x,V) = 3> N (x,V)[X] =: p N (fi,V) 



(4.14) 



(4.15) 



The idea is that <f>^ (x, V ) [A ^] is computable due to the properties of the RPC and the hierarchical structure 
of the cavity field. Let us now disintegrate the Gibbs measure on V x si induced by ( |4.13| l into two Gibbs 
measures acting on V and si separately. To this end, we define the correspondent (random) local free 
energy on V as follows 



¥^\t,x,a,V) :=log 



s N (v) 



exp 



(4.16) 



For a E si, let us define the (random) local Gibbs measure ^N(t,x,a,V) £ ^\(Sn) by specifying its 
density with respect to the a priori distribution as 



4 M) (t,x,a,V) 



(") : =1 Sv(v)Wexp PVnh( m) (u,a) - ^ ) (t ,x,V,a) 



u € Sn- 



Let us define the re-weighting of the RPC by means of the local free energy ( |4. 16| > 

RPC(a) :=RPC(a)exp^v^ M) (f,x,y,a)^) , me/ 

Let us also define the normalisation operation jV : ^#fi n ite O^O ~~ * -^i as 

77(a) 



,(M) 



CC€si, Tj = {r]a)ae.s/ € ^finite 



(4.17) 



(4.18) 



(4.19) 



We introduce the local Gibbs measure ^ N '(t,x,V) £ s$\{V x si) as follows. We equip V x si with the 
product topology between the Borel topology on V and the discrete topology on si. For any measurable 
<H c V x si, let us put 



AM), 



, N >(f,*,V)[*r]:= £ ,yV(RPC)(a)^ \t,x,a,V){veV | (v,a) €<8r}. 



(4.20) 



Let us define the remainder term as 



Sf^ M) (f,jt,V)(d M ,da (1) ) y ^ M) {t,x,V){dv,da {2) ) 
l -(D(2(r-q(a^\a^)))-D(2(r-(u,v) N ))) 



(4.21) 



-D'' M (2(r- ? (a( 1 ),a( 2 ))))( ? (a( 1 ),a( 2 )))-(«,v) w ) 
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Given r £ (0;c/], let us denote 



Define the local remainder term as 



V, 



(r — e;r+e). 



M {M) (r)[x] := lim lim / ^ M, (f,V £ )df, ieJ'(r) 
eJ,+OJVt+~»/o 



(4.22) 
(4.23) 



The main step in the proof of Theorem 2.2 is the following. 



Lemma 4.2 (Comparison with cascades). Given r £ (Q;d], for any x £ <^'(r), as e I +0, and M j" +°° 

^ M \t,x,V e ) = -^ M \r)[x] - j^xt (0 r (M) fe + i) - r (M) (ft)) + 0(e) + ff(\/M), (4.24) 

& {M) (r)[x]>0. (4.25) 
Proof. Fix some r g (0;d]. Using the notation ( |2.21[ ) and smoothness of D on (0; +°°), we have 

VarX( M ) = G r (r) + 0(e), VarA( M , a) = rG' r {r) + 0(e), u£V e , a £ ss/ . (4.26) 

and 

Cov[X(u),X(v)] = G r ((u,v) N ), 

Cov \A(u,aW),A(v,aW)] = G' r (q(a w ,a [2] ))(u,v) N . 

Applying the abstract Gaussian interpolation formula (see, e.g., J5] Proposition 2.5]) to the field Xn and the 
cavity field ( |4.12) , we obtain 



(4.27) 



d B 2 

j- t (p N (t,x,V e (r)) = y E 



y N (t,x,V){&uAa m ) J & N (t,x,V)(dv,da {2) ) 

rVarX(i<)-VarA(H,a)-Cov[X(H),X(v)]+Cov A(u, a {1) ),A(v, a [2) ) ) +0(e). 
Using ( |4.26| ) and ( |4.27| ), we get 

VarX( M ) - VarA(«, a) - Cov [X(w),X(v)] +Cov A(u, a (1) ),A(v, a (2) ) 
= G r (r)-rG' r (r)-(G,.(< ? (a( 1 ),a^))~^(a( 1 ),a( 2 ))G;(< ? (a( 1 ),a( 2 )))) 
- G r ((«,v) w )-G r ( ? (a( 1 ),a( 2 )))-G;.( ? (a( 1 \a( 2 )))((«,v) w - ? (a( 1 \a( 2 )))" 
Comparing ( |2.21| i and ( |2.10| i, we note 

G r (q) - sG' r (q) = D(r) + r (q), q £ R+. 
We have (cf. the proof of |5 , Lemma 5.2]) 

^ M) (f,x,y £ )(d M ,da( 1 ))|^ M) (f,x,y £ )(dv,da( 2 »)(0 r (r)-0,.(^(a( 1 \a( 2 )))) 



(4.28) 



(4.29) 



(4.30) 



= E 



/K(RPC)(da (1 )) / ^(RPC)(da( 2) )(0/ M) (r)-0 r (M) (^(a( 1 ),a (2) ))) 



(4.31) 



= L**(^* + i)-0 r W (<?*)) 
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By ggj) , 



G r ((«,v) iV )-G r ( ? (a( 1 ),a( 2 )))-G;( ? (a( 1 ),a( 2 )))((«,v) i v- ? (a( 1 ),a( 2 ))^ 
= i(D(2(r- ? (a( 1 ),a< 2 ))))-D(2(r-( M ,v) w ))) ( 4 - 32 ) 
-D'(2(r- 9 (c(( 1 U (2) )))(?(a (1 »c( (2) ))-(",v'W- 



Combining ( |4.31| >, ( |4.29| >, ( |4.32| > and ( |4.28| l, we get ( j4.24| l. Due to Remark [O] the function G is convex. 
Therefore, 

G r (( H ,^)-G r ( 9 («( 1 )a (2) ))-G:(?(« (1) ,« (2) ))((« I v) w - ? (« (1) ,«( 2 ))) >0. (4.33) 
Inequality ( |4.25| ) follows from ( |4.33| >. □ 

4.3 Regularisation and localisation 



In this section, we finish the proof of Theorem 2.2 



Lemma 4.3 (Regularisation, well-definiteness). For any x £ ^(r), 

j3 2 " 



lim 



lim lim Mx,V,m - P - J> (e { r M \q k+l )~ 9^ M \q k ) 



x„tl-0 \ £|+u ^ k=l 

Proof. Recall |4TTT) . Given x £ J^'(r), large enough given M > 0, as £ J, +0 and jc„ j" 1 — 0, we have 



(4.34) 



(p^\0,x,V £ ) = ^(M-D'(2(r-q,,)))r + ® N (x,V £ )[A] + 0(e) + 0(1 -x„), 
where A(u,a) := -j= EiLi 5,. (oc)m,-, and {a,} are i.i.d. copies of 



(4.35) 



E 



Using the definition ( |2.10| >, for large enough given M > 0, as x„ 1 1 — 0, we get 

x„(e { r M) ( qk+l )-0^ M) ( qk )) =(M-D'(2(r-q n )))r- l -D(2(r-q n )) + ff(\-x n ). 



(4.36) 



(4.37) 



Combining ( |4.35| l and ( |4.37[ ), we note that the unbounded in M terms in ( |4.34[ ) cancel out and therefore 
g34j holds. □ 



Lemma 4.4 (Localisation, large deviations and cascades). For any x £ 3^l(r), 



Um <pf \o,x,V e ) = inf f/^(0,0)-Ar 

£1+0 AGR l 



(4.38) 



Proof. This is a standard computation (cf., e.g., (U Lemma 6.2]), using the well-known averaging prop- 
erties of the RPC (see, e.g., J5] (5.27)]) and the quenched large deviations principle as is done in ||5] 
Sections 3-5]. □ 



Proof of Theorem \2. 2\ Combining Lemmata |4.2||4.4|and|4.3| we obtain Theorem |2. 2 



□ 



5 Outlook 



Combining the methods of Talagrand ifTBI with Theorem |2.2| we can show that the remainder term in \2. 11) 
vanishes at the saddle-point. This implies that, in fact, the equality holds in ( |2.18[ ). Summarising, we arrive 
at the following result. 



Theorem 5.1 (Parisi-type formula). In the case of the product state space \2.5\, for all j3 € R 



p(j5) = sup inf &(l3,r,g)[x\, almost surely. 



(5.1) 
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Parallel to the product state space ( |2.5| >, one can consider the rotationally invariant state space: 

S N :={ueR N :\\u\\ 2 <lVn}, L > 0. (5.2) 
In this case, we assume that the a priori measure jj.^ £ ^#fi n j te (Sn) has the density 



d;U 
dl 



/V 



(u):=«p , u=( Ui )f =1 eR N , f:R^R (5.3) 



vi=l 



with respect to the Lebesgue measure A on R N . Let the function / be of the form f(u) := /ijm — ti2U 2 , 
where hi £ R and h% £ R+ are given constants. Let us note that in case ( |5.2) , c/ = L 2 . 

In the case of the rotationally invariant state space ( |5.2) , one can obtain a more explicit representation 
for the Parisi functional d2J3), which does not require any regularisation. Given x £ SE(f), define q msK := 
q max (x) := supjg £ [Q;r] : x(q) < l}. Consider the Crisanti-Sommers type functional (cf. (5J (A2.4)] and 

as (47)D 



\og(r-q max ) + / +/t] / *(0)d$- 

Jq w J (5.4) 

6 2 / /•'?max \ 

T D'(2(r- 9 ))jc(?)d 9 , iel(r). 



o 

By reducing the case of the rotationally invariant state space to the product state space case using a large 
deviations argument (an idea exploited in |[T5l ). one arrives at the following result. 

Theorem 5.2 (Fyodorov-Sommers formula). In the case of the rotationally invariant state space \5.2\ , for 
all j3 £ R+, h\ £R, li2 £ R+, there exists unique r* £ [0;d] and unique x* £ 3£{f) such that 

p(/3)=max min f ^(j3,r)[x] = ^Q3,r* )[**], almost surely. (5.5) 

re[0;d]xe3r{r) 



The proofs of Theorems |5.1| and |5.2| are beyond the scope of this short communication and will be 
reported on elsewhere. 



Remark 5.1. The Crisanti-Sommers type functional ( |5.4| i corresponds to the a priori distribution ( |5.3| l, 
which represents the linear combination of linear and quadratic external fields. Formula [8, (47)] was 
derived under the assumption of the quadratic external field, whereas formula [6. (A2.4)] was obtained for 
the spherical SK model with the linear external field. 



Remark 5.2. The explicit form of the functional ( 5.4 1 assures that it is strictly convex with respect to 
x £ S£(f). In contrast, convexity of the functional ( |2.13 I is (to the author's best knowledge) open, see HI2V 
and Theorem 6.4] for partial results. 

A Characterisation of the correlators 

We recall some facts about high-dimensional Gaussian processes with isotropic increments. The following 
result can be found in the work [ 19| of A.M. Yaglom (see also |20|). 

Theorem A.l. IfX is a Gaussian random field with isotropic increments that satisfies ( |2.1| l, then one of 
the following two cases holds: 

1. Isotropic field. There exists the correlation function B : R + — > R such that 

E[X N (u)X N (v)]=B(^\\u-v\\ 2 1 j, M ,veE/v, (A.l) 
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where the function B has the representation 

B(r)=c + jJ exp(-r 2 r) v(df), (A.2) 

where cq £ R + is a constant and V € .^finite is a non-negative finite measure. In this case, the 
function D in ( |2.1| l is expressed in terms of the correlation function B as 

D(r)= 2(5(0) -B{r)). (A.3) 

2. Non-isotropic field with isotropic increments. The function D in \2.\\ has the following representa- 
tion 

D(r) = J^ [l-exp(-r 2 r)]v(dr)+A-r, reR+, (A.4) 

where A £ R + is a constant and V £ ^#((0; +°°)) is a O-finite measure with 

r+°° t 2 v(dt) 



t 2 + l 



< oo. (A.5) 



Remark A.l. In Theorem A.l assuming co = 0, ca.se[7]i.s sometimes referred to as the short-range one 
which reflects the decay of correlations: B(r) J, +0, as r +°° This fact follows from the representa- 
tion ( |A.2| >. Correspondingly, case [2] is called the long-range one, since here, assuming X(0) = 0, the 
correlation structure is 

E[X N (u)X N {v)} = - (Dn(\\u\\1)+D n (\\v\\1)-D n (\\u-v\\1)) , u,veR N . (A.6) 



Equation ( |A.6| > in combination with the representation ( |A.4| > implies that the correlations of the field Xn do 
not decay, as \\u — v| — > +°°. 

Remark A.2. Theorem \A.1\ implies that the function D appearing in ( |2.1| l is necessarily concave, infinitely 
differentiable, and non-decreasing on (0;+°°). 
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